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Abstract 

The principal theme of this paper is that anomalously slow, super- Arrhenius relaxations in glassy 
materials may be activated processes involving chains of molecular displacements. As pointed out 
in a preceding paper with A. Lemaitre, the entropy of critically long excitation chains can enable 
them to grow without bound, thus activating stable thermal fluctuations in the local density or 
molecular coordination of the material. I argue here that the intrinsic molecular-scale disorder in 
a glass plays an essential role in determining the activation rate for such chains, and show that 
a simple disorder-related correction to the earlier theory recovers the Vogel-Fulcher law in three 
dimensions. A key feature of this theory is that the spatial extent of critically long excitation 
chains diverges at the Vogel-Fulcher temperature. I speculate that this diverging length scale 
implies that, as the temperature decreases, increasingly large regions of the system become frozen 
and do not contribute to the configurational entropy, and thus ergodicity is partially broken in the 
super- Arrhenius region above the Kauzmann temperature Tr. This partially broken ergodicity 
seems to explain the vanishing entropy at Tx and other observed relations between dynamics and 
thermodynamics at the glass transition. 
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I. INTRODUCTION 

In a preceding paper Lemaitre and I explored the hypothesis that anomalously slow, 
super- Arrhenius, relaxation rates in glassy solids can be understood by assuming that transi- 
tions between the inherent states of such materials are enabled by thermally activated chains 
of small molecular displacements. (See |2j for a summary of research in a wide range of topics 
related to the dynamics of glassy materials.) More specifically, we developed a model of the 
spontaneous formation of shear-transformation-zones (STZ's) by thermal fluctuations in the 
absence of driving forces. STZ's are localized irregularities in the density of molecules, or in 
near-neighbor molecular correlations, that undergo irreversible rearrangements during shear 
deformation. P, U, 0] We visualized their formation as the glassy analog of the formation of 
vacancy-interstitial pairs. At low temperatures, such pairs must become well separated in 
order to be thermodynamically stable against recombination; and the state that includes 
such a well separated pair is inherently distinct from the state in which the pair is absent. 
In our model, the activation energy for forming a stable pair is the free energy - including 
the entropy - of a chain of molecular displacements that moves the "interstial" away from 
the "vacancy" and is just long enough to be marginally unstable against further growth. In 
short, we solved a nucleation problem in which the relevant reaction coordinate is the length 
of the chain. The entropy associated with different chain configurations is a measure of the 
number of routes across the activation barrier (see [6(), and therefore plays a central role in 
determining the transition rate. 

The ana^is presented in Q wa 8 at best o„. y pa rtly successful. The excited chain, 
seemed to have qualitatively the right properties to predict a diverging viscosity at a non-zero 
temperature To. In order to produce a Vogel-Fulcher law, however, we had to assume that the 
chains were restricted to lie on two-dimensional surfaces, perhaps the interfaces beween the 
frustration-limited domains of Kivelson et al [3], or boundaries within the mosaic structures 
proposed by Kirkpatrick, Thirumalai, Wolynes and others. P. lil llfl| This picture may indeed 
deserve further investigation; I shall return to it briefly in discussing the thermodynamics of 
these systems. The main purpose of the present paper, however, is to argue that the missing 
feature of the model discussed in [l| may be the spatially disordered environment that is 
intrinsic to any glassy material. 

The problem of computing the activation energy of excitation chains is related - but not 
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exactly equivalent - to the problem of computing the statistics of self-avoiding random walks. 
That problem, in turn, is approximately equivalent to solving a diffusion equation in a self- 
consistent repulsive potential. In [lj, we approximated the solutions of those related prob- 
lems by adapting Flory's method for calculating excluded- volume effects in polymers. fill] 
I do essentially the same thing here, but solve a (non-conserving) diffusion equation in a 
random potential whose spatial disorder is that of the configurational degrees of freedom 
of the glass. It is well known that diffusion is constrained in disordered systems; thus the 
statistics of excitation chains must be determined by the interplay between two competing 
effects: swelling of the excitation region due to excluded volume, and contraction due to 
disorder. The result of this competition is the Vogel-Fulcher law in three dimensions. 

Before entering into the details of this analysis, a few remarks about how it is related 
to previous work in this field are in order. The scientific literature abounds with attempts 
to explain super-Arrhenius behaviors in glasses. Some of the most important and influen- 
tial of these papers, notably that of Adam and Gibbs Q], are almost half a century old. 
More recent work along these lines has been carried out by Wolynes and his collaborators 
aSQand by o th e r invests a s renewed in Q. (See Q for a perceptive analysis 
of the 'Adam-Gibbs-Kirkpatrick-Thirumalai- Wolynes scenario.") All of these authors make 
one fundamental assumption with which I concur - that the super-Arrhenius relaxation pro- 
cesses are intrinsically nonequilibrium phenomena. By this I mean that the super-Arrhenius 
formula pertains s rictly to the rate at which a glassy material makes transitions between 



its inherent states Jjj |JJ , 



13], and not just to the statistics of those states themselves. This 



assumption sets us apart, for example, from Cohen and Grest who, in their classic paper 
[3| > attributed super-Arrhenius behavior to percolation of liquid-like regions in equilibrated 
states. An especially relevant counterexample is the long-standing convention in engineering- 
oriented papers on metallic glasses and other amorphous systems, where it is assumed that 
the equilibrium density of flow defects is a super-Arrhenius function of temperature. For 



example, see 

l2d, |2l|, |2^, Q, Q • 1 sha11 not dwe11 on this 

issue, but it needs to be taken 
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into account when reading that literature. 

The principal way in which the super-Arrhenius processes are relevant to thermodynamic 
equilibrium is that they determine the validity of the ergodic hypothesis. It is in this 
connection that I disagree with the approach of Wolynes et al, which is based on the idea 
that the activated transition states consist of droplets of an entropically favored liquid phase. 
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Such a phase seems to be precluded by the same, extremely slow relaxation rates that all of 
these authors are trying to compute; it cannot exist on experimentally relevant time scales, 
especially not as an ephemeral transition state. Rather, the activated chains discussed in 
and in the present paper might provide a better model of the transition mechanism because 
they are ordinary thermal fluctuations occurring within single inherent states. Moreover, 
the excitation-chain model involves only short-ranged interactions between the constituent 
molecules; it draws none of its ingredients from infinite-range mean-field models or the like. 
Thus, at least in principle, it should be possible to use it to develop microscopic descriptions 
of various kinds of glassy systems. 

The scheme of this paper is as follows. Section II contains a summary of basic ideas 
and principal predictions of the excitation-chain dynamics. Section III is devoted to the 
mathematics of the disorder problem in glassy systems and the derivation of a formula used 
in Section II. Readers who are not interested in such details may skip that part of the pa- 
per; but the analysis described there addresses some technical issues that distinguish this 
nonequilibrium situation from otherwise similar problems. In Section IV, I propose an in- 
terpretation of the preceding results that seems to explain the striking relationships that 
have been found to exist between the dynamics and thermodynamics of the glass transition 
- that is, the apparent equivalence of the Vogel-Fulcher and Kauzmann temperatures, the 
success of the Adam-Gibbs theory, and the approximate proportionality between the dy- 
namic fragility and the jump in the specific heat at the glass temperature. I close with a 
list of unanswered questions. 



II. EXCITATION-CHAIN DYNAMICS 



All of the discussion that follows is based on a picture of a glass as a supercooled liquid in 
which configurational rearrangements have become very much slower than the thermal vi- 
brational motions of the molecules within their local environments, i.e. within their "cages." 
The goal of the excitation-chain theory is to compute the rates of configurational rearrange- 
ments, and to do this using a mechanism that involves just the rapid thermal fluctuations 
that ultimately must drive those motions. An essential element of this picture is the notion 
of frustration - that the energetically most preferred local configurations of the molecules do 
not fit together to fill space, and therefore that there is a geometrically necessary population 
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of somewhat higher energy "defects" where the local coordinations are not the most pre- 
ferred ones. We may think of the distribution of these relatively populous, necessary defects 
as merging, at yet higher energies, with a much smaller population of shear-transformation 
zones (STZ's) and other localized irregularities in the density and/or molecular arrange- 
ments. The anomalously loose STZ-like defects apparently govern the response of the system 
to external driving forces. They are the source of the observed super-Arrhenius behaviors, 
and therefore are of special interest here. 

The relevance of one-dimensional excitations to the nonequilibrium behavior of glassy 
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s is supported by the molecular dynamics simulations of Glotzer and colleagues, 



271 ] which showed that transitions between inherent states in glass-forming liquids 



take place via motions of stringlike groups of molecules. The well documented existence 
of force chains in granular materials Q seems to be further evidence in favor of the idea 
that forces and displacements are transmitted primarily along one dimensional structures in 
noncrystalline systems. 

As in p|, consider just the spontaneous creation of an STZ, that is, its creation due to 
a thermal fluctuation in the absence of external driving; and assume that such an event is 
roughly similar to the formation of a "vacancy-interstitial" pair, followed by displacements 
of atoms along an excitation chain. Lemaitre and I visualized this chain as a thermal fluctu- 
ation in which a linear array of momentarily loosened atoms undergoes small displacements, 
effectively moving the "vacancy" and the "interstitial" far enough away from each other 
that they do not quickly recombine. More precisely, the transition state for this activated 
process is a momentary thermal excitation of the system that enables a chain of molecular 
displacements just long enough that it is as likely to grow as to decay. No molecule has 
yet moved fully out of its cage in this transition state. Excitation chains smaller than the 
critical size, like subcritical liquidlike clusters in a supercooled vapor, with high probability 
just disappear, leaving the system unchanged. Once the chain exceeds this critical size, 
however, the vacancy and the interstitial become uncorrelated with each other, and each 
finds its own stable position in a new inherent state. The time taken by such a transition, 
once it occurs, is roughly the time during which the chain fluctuates between different near- 
critical lengths and configurations as it passes across the activation barrier. This time may 
be very long compared to an oscillation period for a molecule in its cage, but it is very short 
- essentially instantaneous - in comparison to the inverse of the super-Arrhenius rate at 
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which these transitions occur. In other words, at temperatures near the glass transition, 
excitation chains are very rare and relatively brief events. 

The challenge is to compute the probability per unit time for formation of an excitation 
chain. Consider a chain of length N, measured in units of a characteristic molecular spacing 
which, without loss of generality, can be set to unity; and suppose that the chain occupies a 
roughly spherical region of radius R in a three-dimensional space. Throughout the following, 
R is the expected distance from the origin to the last, iV'th link in the chain; but it seems 
reasonable to assume that the radius of the occupied region is the same as R up to an 
unimportant geometrical constant. 

The required formation probability for the pair plus the excitation chain is the product of 
a Boltzmann factor containing the activation energy, multiplied by the number of configura- 
tions of the chain with length N and extension R. It is conventional to write the logarithm 
of this probability in the form — AG(N, R)/ks T, where AG(N,R) is called the activation 
/ree-energy because it includes something like an entropy. It will become clear in Section III 
that this quantity is not exactly a free energy in the conventional sense. 

AG(N, R) consists of several parts: 

AG(N, R) = AGoo + Ne + E int (N, R)-k B T In W(N, R). (2.1) 

The first term, AGqo, is the bare activation energy, that is, the energy required to form the 
"vacancy" and the "interstitial." AGoo becomes the ordinary Arrhenius activation energy at 
high enough temperatures, say T > T4, where N must vanish because the chain is no longer 
needed to stabilize the excitation. (See remarks at the end of this Section.) In equilibrium 
situations even below Ta, AGoo is the activation energy that occurs in a Boltzmann factor 
for determining the population of these defects; and detailed balance requires that this 
Boltzmann factor be the ratio of their creation and annihilation rates. Different processes, 
associated with different kinds of defects, will have different values of AG^. On the other 
hand, in nonequilibrium situations below T4, the chain dynamics near the glass transition 
may be much the same for different defect-related mechanisms because the chains involve 
large numbers of molecules. Thus the super- Arrhenius behavior may be common to a variety 
of different relaxation phenomena in a single material. 

The remaining terms in Eq. (j2.1j) describe the excess free energy of the chain. The second 
term on the right-hand side of Eq. (j2.1j) . iVeo, is the bare activation energy of the N links of 
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the chain, unmodified by entropy or the self-exclusion effect. The average energy per link, 
eo, is a measure of the elastic stiffness of the molecular environments, that is, the energy 
required to move two molecules far enough apart from one another to allow a third to pass 
between them. 

The third term, E int , makes it energetically unfavorable for the links of the chain to lie 
near one another. In principle, this exclusion effect might be included directly in a sum over 
self-avoiding random walks, perhaps using a nonperturbative method like that described 
by Edwards. j3] For present purposes, however, it seems better not to be so ambitious, 
especially since the exclusion forces in this case are likely to be long-ranged. The thermal 
fluctuation that loosens the molecules along the chain must push molecules closer together 
at points away from the chain, thus producing an extended repulsion. As a result, Flory's 
mean-field approximation may be more accurate here than it is for polymers. As in the 
Flory interaction energy is proportional to the square of the string density multiplied by the 
volume occupied by the string. That is, 

N 2 

E int (N,R) = k B T mt —, (2.2) 

where dimensionless geometric factors have been absorbed into the definition of T int . Note 
that this approximation makes sense only in the limit of large N. The exclusion effect must 
disappear for short chains - an important complication that will be discussed later. 

In the last term on the right-hand side of Eq. (j2.1j) . W(N, R) is a sum over chain config- 
urations. Evaluating W(N,R) is the crux of the present analysis. In lj, we wrote 

\nW(N,R)&vN-^, (2.3) 

where exp (u) is the number of choices that a walk can make at each step and exp (—R 2 /2N), 
the free-diffusion factor (up to a normalization constant), is the a priori probability for a 
chain of length N to occupy a region of radius R. We then minimized AG(N, R) with respect 
to R and maximized it with repect to N. That is, we found a saddle point of AG(N, R), 
which we identified as the activation energy. In three dimensions, our result was 
AG*(T) AG(N*,R*) AG^ T 3 ^T^ 

T"^ = \ « -jT 1 + C ° nStant X IT T 2 ' 4 

k B k B k B (T-Tq) 1 ^ 

where T = eo/{i>k B ), and (R*, N*) is the location of the saddle point: 

,rp x 1/5 7^3/4 7^1/2 



^oc(^) (AT*) 3 / 5 ; APoc ( ^_^ /4 . (2.5) 



The exponent 1/4 in Eq. (j2.4j) was clearly too far from the Vogel-Fulcher formula to be 
consistent with experiment, and thus we looked at its analog in two dimensions where that 
exponent turns out to be unity. (Dimensionality entered only via the interaction term.) 

The thesis here is that the missing ingredient of the preceding analysis was the intrinsically 
disordered, glassy environment in which the excitation chains occur. This disorder appears 
in the activation energies for individual links of the chain, and the variation of these energies 
reflects the structural disorder in the underlying molecular configurations. To incorporate 
this disorder into the evaluation of AG(N, R), write the activation energy at the position of 
the n'th link of the chain, say r n , in the form eo + ksT ip(r n ), where ip(r n ) is a dimensionless 
random variable. If the disorder is uncorrelated from site to site, and eo has been chosen to 
be the average activation energy per link, then 

< (p(r n ) >= 0; < y?(r n ) <p(r m ) >= j(T) 5 n>m , (2.6) 

where the angular brackets denote a statistical average over realizations of the disorder. The 
function j(T) is the strength of the disorder associated with the geometrically necessary, 
frustration-induced, configurational defects discussed at the beginning of this Section. These 
defects must pervade the system, and their density should not be a strongly varying function 
of temperature. Let this defect density be n d per molecule, and further suppose that the 
variations in eo are of the order of eo = v ks Tq. Then, 

7(T)= 7 o(^) 2 ; lo = v 2 n d . (2.7) 

Equation (j2.7J) has several important implications. Note that 7(T) is a decreasing function 
of the bath temperature T. This feature reflects the fact that, at lower temperatures, 
the chains are more tightly constrained to lie in the minima of the random field <p(r n ), 
and thus the effective coupling between the disorder and the chains is greater. Note also 
that 7 may be of order unity, especially in the neighborhood of T = T . The density of 
geometrically necessary defects is likely to be substantial and, at lower temperatures, the 
fact that the activation energy changes rapidly from site to site may invalidate the simple 
approximation made in Eq. (J2.7|) . If 7 is large, then a weak-coupling expansion of the kind 
described in Section III would not be quantitatively accurate. However, since this theory is 
meant primarily to be a plausibility argument in favor of the disorder hypothesis, it seems 
reasonable to assume that it gives at least a qualitatively correct answer. 
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The crucial prediction of the disorder theory developed in Section III is that the free 
diffusion factor, exp (—R 2 /2 N) in W(N, R), becomes negligable in comparison with a local- 
ization factor exp (—7i r yR/2) for large N and R. Eq. (j2.3j) becomes 

In W(N, R)^uN- n 1 ^ R. (2.8) 



Inserting Eq. (|2.8j) into Eq. (|2.1|l and computing the saddle point, I find: 

AG*(T) = AG(N*,R*) ^ AGoo /vr\ 3 / 2 jiXf^T^T^ 

k B k B k B \6J u{T-T ) ' 1 ' ' 

where 

/7T\ 3 / 2 7 m 3 / 2 T 3 / 2 T 1 ( 2 

N " Ki (d J(t-^ " ' (2 - 10) 

and 

ir-.-( 6r -' V V)'/'~ 2 (lY* (211 ) 

To justify neglecting the diffusion factor, note that the first expression for R* in Eq. ()2.11|) 
implies that R* 2 /N* is of order unity, while 7 R* grows like N* 1 ^ 2 for large iV*. Thus, the 
disorder effect restores the Vogel-Fulcher result in three dimensions, and restores agreement 
with experiment near T without invoking frustration-limited domains or mosaic structures. 

These results, as they stand, do not account for the transition between liquidlike and 
solidlike glassy behavior that occurs when the excitation chains disappear at T = Ta- The 
fact that Ta is a well defined temperature is supported both by experimental evidence such 
as that shown, for example, in Fig.l of 2^, and by the excitation-chain idea itself. At 
temperatures below Ta, a molecule that makes a thermally activated jump to a neighboring, 
energetically unfavorable "interstitial" position is most likely to jump back to its original 
position in its next thermally activated transition. At higher temperatures, on the other 
hand, the Boltzmann probability that favors recombination can be compensated by an en- 
tropic factor of the form exp(ui) iv\ = u), which counts the number of allowed jumps that 
move the interstitial further away from the vacancy. In other words, to evaluate Ta, we 
must compute the temperature at which the critical length of a chain is A^* = 1. Just as the 
exclusion energy and the disorder are relevant when A^* is large, the local environment of the 
vacancy and interstitial must be relevant in computing Ta- An accurate evaluation of E int 
for small and R, or preferably a non- mean-field theory of the small-chain limit, will be 
needed in order to construct a quantitative theory of the transition between super- Arrhenius 
and Arrhenius behaviors. That problem is beyond the scope of this paper. 



III. PERTURBATION THEORETIC ANALYSIS OF THE DISORDER EFFECT 



rr(N)=R 

g Q {N,R,0) = / 8r(n) exp 

Jr(o)=o 









(£) dn 







(3.3) 



This Section contains a somewhat old-fashioned demonstration of how the exponential 
decay law in Eq. (J2.8|) emerges in a self-consistent perturbation-theoretic approximation. 
This approximation is srictly valid only in the limit of small 7. A better calculation, using 
more advanced techniques, should be feasible. 

The function W(N,R) in Eq. (j2.1j) is a sum over chain configurations. In [lj], Lemaitre 
and I approximated it by writing 

W(N, R) « W (N, R) = e Nu Q (N, R, 0), (3.1) 

where Go(N, R, 0) is the three dimensional diffusion kernel with a source at R = 0: 

^o(iV,i?,0) = ^ F ,e-^^. (3.2) 

It is useful to think of this kernel as a Wiener integral over continuous paths r(n), where n 
is a continuous variable running from to N, and r(0) = 0, r(N) = R. That is, 

rr(N)=R 
/r(0)=0 

(See, for example, the classic review article by Gel'fand and Yaglom 30].) In this functional 
form, it is clear that the sum over paths is a version of the desired sum over chain configu- 
rations. The Gaussian exponential factor, i.e. the Wiener measure, constrains the paths to 
be connected, one-dimensional objects embedded in a three-dimensional space. This contin- 
uum approximation is convenient analytically and is perfectly accurate so long as it does not 
make much difference that the steps along a chain have length unity (the molecular spacing) 
and that dn = 1. 

One problem here is that, strictly speaking, the glassy disorder is not consistent with the 
continuum limit. It requires that each unit step along a chain, say the j'th, have its own 
extra weight factor exp (—ipj) in the sum over configurations. In the face of this difficulty, 
it is easiest to solve a slightly different problem. Let the chains remain as continuous one- 
dimensional objects, with no explicit links of finite length; but suppose that they exist in the 
presence of a random potential (p(v), defined over the whole volume of the system spanned 
by the continuous variable r; and suppose that, in analogy to Eq. (j2.6|) . 

<<^(r)>=0; <(p(r)(p(r')>=j5(r-r'). (3.4) 
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The assumption of delta-function (white noise) correlations for the random field </?(r), implies 
that the molecular length scale (unity) on which this field varies is very much smaller than 
any other relevant length scale in the problem, especially R. Also, assume that the delta- 
correlations arise from a Gaussian distribution over the <p(r). That is, the average of any 
function A({ip}) is 

(A) = f 6<p(r) exp -^J V 2 (r)dr A({<p}), (3.5) 

where the functional differential S(p(r) is defined to include a normalization factor so that 
<!>=!. With these assumptions, replace Qq in Eq. ljH.ljl by 



g(N,R,0,{<p}) 



r(N)=R 
r(0)=0 



Sr(n) exp 



N 



dn 



(3.6) 



The weight function W(N,R) must then be obtained by computing (S(N, R, 0, {<£>})) ac- 
cording to Eq. lj3.5jl . 

This model is now in familiar territory. It is closely related to the quantum model for 
a single particle moving in a random potential. G{N, R, 0, {f}) satisfies the differential 
equation 

^-^ + y?(R) g(N,K,0,{ V }) = 5(K)5(N), (3.7) 

Thus Q is the Green's function for the imaginary-time (N) Schroedinger equation for a 
particle moving in in a delta-correlated random potential y?(R) of strength 7. 

An apparently more closely related problem is that of polymer chains in disordered media. 
In this case, the competition between strong disorder and strong self-exclusion has been 
e Xplore d, for e Xampl e, . Q Q . LanaUre lad I US ed the poller m a,o gy in our earlier 
paper jl|. This analogy does deserve to be explored further; but it may not be quite so 
close as it appears. In the polymer problem, one must compute the disorder average of 
the thermodynamic free energy in a quenched system. In the present case, the disorder is 
also quenched in the sense that it is fixed and independent of the chain configuration; but 
the chain is not a pre-existing entity. The goal is to compute the disorder average of the 
probability that the chain appears during thermal fluctuations in the system; and thus the 
required quantity, (Q(N,R,0, {<£>})), more nearly corresponds to an annealed average. 

The next step in this analysis is to construct a perturbation expansion for Q. To do this, 
it is easiest to work with the Laplace transform of Q: 

~ poo 

Q(w,R,0,{<p})= dNe- wN g(N,R,0,{<p}), (3.8) 
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which satisfies 

w ~ + p(R)] R, 0, M) = S(R). (3.9) 

It also is necessary to renormalize w by writing w = wo + w' and using wq + y> as the 
perturbation. The renormalization constant il>o will be chosen as in old-fashioned particle 
field theory to cancel a formally (but in this case not really) divergent integral. 

Now expand Q in powers of this perturbation and average the expansion term by term over 
ip. (The basic idea for averaging a perturbation expansion over a random potential goes back 
to Kohn and Luttinger [33. That procedure, and most of the other techniques used in the 
next paragraphs, are described in Chapter 4 of Mahan's book on many-particle physics. [34]) 
Averaging over the disorder restores translational symmetry, which makes it natural to work 
with Fourier transformed functions. Let Q(w, k) denote the Fourier transform of the average 
of G(w, R, {ip}). The calculation of Q(w, k) then follows completely conventional lines. The 
re-summed perturbation expansion for Q(w, k) has the form 

S«k)= fc2 / (3.10) 

where here the self-energy E is a function only of w'. In general, E would be fc-dependent; 
but that dependence would arise from short-range spatial correlations in y?(r), which are 
assumed to be absent. In fact, the molecular length scale, taken here to be unity, provides 
a necessary short wavelength cutoff. In Eq. ()3.1()j) . the relevant values of k are of order 
1/R <C 1, so it is appropriate (again in the limit of large R and N) to make a small- A; 
approximation - except, of course, when one encounters short-wavelength divergences. 
The lowest order, self-consistent approximation for E is: 

r rl 3 k' ~ 

E(w')^w Q - 1 j-^- 3 g(w',k'), (3.11) 

Note that it is Q and not the unperturbed propagator that appears on the right-hand side 
of Eq. (j3.11|) . This self-consistent approximation ensures that E(w') and Q(w',k) have the 
same analytic structure in the complex w' plane, a condition that is known to be true on 
very general grounds and which turns out to be essential for present purposes. 
Inserting Eq. (j3.1U|) into Eq. (j3.11|) yields 

d k ~ , / . /x kmnx 1 



Q[w ^ k') = ^- = Jw> + EK). (3.12) 
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Here, k max is the anticipated short-wavelength cutoff that is needed in order to evaluate an 
w'- independent divergent integral over k'. The renormalization constant conveniently can 
be chosen to cancel this nominal divergence: 



k 

""max i n -, o \ 

to = 7— j-, (3.13) 



7T 

so that Eq. (|3.11|) becomes 



7 

Z.\W) = - 

n 

and therefore 



Z(w') = -^=^w> + Z(w>), (3.14) 



TPS 2 



'l + ^T^ + 1 
7T 7 Z 



4 

A second shift in the Laplace variable, w' = —n 2 7 2 /8 + w" leads to 



(3.15) 



2 



SK, fe) = - : !_ -2 , (3.16) 



The last step in this part of the mathematical development is to invert the Fourier and 
Laplace transforms. 

dw" f+™ d 3 k e ikR+Nw" 



(G(N, R)) = e~ NAeo ^ kBT 



2m J-ao (2tt)3 £ + ^ + ^ 2 



- (2 rf) 3/ 2 v 1 + I75j( 1 + ^ e "™ _1 ' (3 ' 17) 

The final result is accurate up to corrections of relative order 1/N. All of the shifts of the 
Laplace variable (renormalization corrections) are combined here into a shift of the bare 
excitation energy e : 

Ae (T) , . k max vr 2 7(T) 2 , 

= _ 7(T) ^£ + n ) (318) 

This expression is apparently the beginning of a series in powers of 7 that may - or may 
not - be convergent. The main effect of Aeo(T) is to shift T ; and, because we have no a 
priori estimate of the unshifted To, it is easiest to omit this term altogether. On the other 
hand, if 7(T) is large, its temperature dependence will be important in fitting experimental 
data near To. In what follows, I omit Aeo(T), and note simply that Eq. (|3.17j) exhibits the 
expected behavior specified in Section II, specifically, the factor exp {—n^R/2). 

The various shortcomings of this perturbation-theoretic result point to needs for further 
investigation. The approximation used here is unlikely to be accurate enough for exploring 
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the excitation-chain model quantitatively if the values of 7 are as large as anticipated. 
Also, the calculation properly should include the exclusion effect along with the disorder, 
instead of dealing with each of these effects separately as if they were decoupled from one 
another. And an accurate theory should account more carefully for the molecular-scale 
features that must play an important role near Ta, where the chains become short. What 
theoretical techniques might be effective for solving the large-7 versions of this model? 
Perhaps diagrammatic techniques can be helpful in computing systematic corrections to the 
localization factor exp (— ir 7 R/2). But perturbation expansions are notoriously incapable of 
predicting some of the most interesting behaviors in systems of this kind. The only possibly 
relevant, nonperturbative, strong-coupling approach with which I am familiar is my 1966 
calculation, in collaboration with Zittartz, of the density of states in an electronic impurity- 
mnd tail, where the result has an essential singularity as a function of the disorder strength. 
35I ] That calculation, however, dealt specifically with the statistics of deeply bound states 
in a white-noise potential, and not the diffusive motion of a particle in that potential. I 
have not yet found a way to apply that nonperturbative method to the present situation. 



(See 



El, 



2j for some ideas about non-perturbative calculations in the related context of 



polymer chains in disordered media. 



IV. THERMODYNAMICS 



The excitation-chain theory makes no assumptions about the underlying thermodynamic 



properties of the glassy materials in which the excitations occur. In this sense, it is c 



36 



37 



38, 



oser 



to the kinetically constrained models of Fredrickson, Andersen, and others 
than it is, for example, to Derrida's random-energy model 40] or Wolynes' hypothesis of 
a random first-order phase transition. j|] The kinetically constrained models have dynamic 
properties that look very much like those of glasses - i.e. dramatic slowing of relaxation 
rates similar to that found in the present theory. One might even postulate that STZ-like 
defects, containing excess free volume, could serve as the facilitating sites introduced in 
EJ. (See also the work of Garrahan, Chandler, and coworkers 0, 42, I^J for a related 
point of view about the roles of facilitating sites.) Like the present model, the kinetically 
constrained models require no thermodynamically singular behavior. It might be interesting 
to learn whether they would exhibit anything like a Kauzmann phenomenon if they were 
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used in numerical simulations of differential scanning calorimetry experiments. 

On the experimental side, decades of careful thermodynamic measurements indicate that 
the configurational entropy of glassy materials drops linearly toward zero at the Kauzmann 
temperature Tk - a behavior that is known to be consistent with equilibrium statistical 
mechanics only for model systems with long-range and usually built-in random interactions 
between their constituent elements. There are striking relations between the Kauzmann 
thermodynamic phenomenon and the dynamics of real glasses. The Kauzmann temperature 
Tk seems to be very close to, and possibly exactly the same as, the Vogel-Fulcher tem- 
perature To. The Adam-Gibbs theory(l2j], which says that the logarithm of the viscosity 
is proportional to the inverse of the configurational entropy near T = Tk, seems to be 
consistent with a significant body of experimental data. And there are indications of at 
least a correlation, and perhaps a direct proportionality between the measured jump in the 
configurational specific heat at the glass transition and the dynamic fragility in a wide range 
of glassy materials. [9| I find myself deeply puzzled by this situation. It is hard for me to 
believe that long-range, mean-field models can accurately predict the dynamical behavior 
of glasses consisting of small molecules with short-range interactions. Moreover, as men- 
tioned in the Introduction, I am not convinced that the connection between dynamics and 
thermodynamics has properly been established for mean-field models that exhibit a ther- 
modynamic Kauzmann phenomenon. On the other hand, the thermodynamic observations 
must be taken very seriously. 

The excitation-chain theory provides one clue that may point toward a resolution of the 
dilemma. This purely dynamic mechanism predicts a length scale, R*(T) in Eq. (|2.11j) . that 
diverges like (T — T ) _1 near T and vanishes at T4. An isolated region of the system that 
is smaller than R* would not be large enough to sustain critically long excitation chains; 
therefore its constituents would be frozen and unable to make thermally activated transi- 
tions between different inherent states. (The molecules in such a region, of course, would 
continue to undergo rapid thermal fluctuations within their cages; but these fluctuations 
would be unable to activate configurational rearrangements.) Regions larger than R*, how- 
ever, would support critical excitations and be able to change their sizes and shapes via 
molecular rearrangements, albeit very slowly at low temperatures. If such regions exist in 
any meaningful way, either as observably bounded domains or simply as slowly fluctuating 
volumes in which the correlations are extremely long-lasting, then their characteristic sizes 
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would be proportional to R*. Smaller regions would be frozen and would be able to grow 
only by coming into contact with larger, unfrozen regions; the latter regions would shrink 
to increase the entropy of the system whole. 

The question of whether spatial heterogeneities occur in glass forming materials seems not 
yet to be clearly resolved. See reviews by Sillescu |45[ and Ediger |46j|, who argue in favor of 
heterogeneity. In a recent paper, Berthier et al j^J point out that heterogeneity in glasses can 
be detected by measuring multi-point dynamic susceptibilities, and show experimental data 
indicating that heterogeneities exist. Shi and Falk 48|, in STZ-related molecular dynamics 
simulations, have seen spatial patterns apparently similar to what I propose here. They 
find that their two dimensional Lennard- Jones glass, when annealed, consists primarily of 
domains in which the molecules are strongly correlated in low-energy configurations, and 
that these strongly correlated domains are separated by interfaces containing higher energy 
defects. 

To solve the thermodynamic puzzle, I postulate that the only unfrozen configurational 
degrees of freedom of this system exist on the boundaries of strongly correlated regions of 
size R*. In other words, I conjecture that the configurational entropy per molecule of the 
system as a whole is proportional to the surface-to- volume ratio, 1/R*, of such regions. 
Specifically, 

s(T)~^= 1 -( 6 -) 1/2 ^ T ~ T °) (4 1) 

where s c (T) is the average configurational entropy per molecule in units of and So is the 
configurational entropy per unfrozen molecule. The entropy Sq is a measure of how free the 
boundary molecules are to rearrange their positions, orientations and, in the case of complex 
molecules, their internal configurations - freedoms that presumably are lacking inside the 
frozen regions. 

Here, and in what follows, I assume that T is sufficiently close to T that the long-chain 
approximation is accurate, and that the super- Arrhenius part of Eq.(j2HJ) dominates the 
Arrhenius part, AG^. Equation (J4.1j) should not be taken literally out to temperatures in 
the neighborhood of T4, where - as in the case of AG*(T) - some correction for vanishingly 
short chains will be needed, and where s c (T) must join smoothly to the configurational 
entropy of the liquidlike state. 

Equation (J4.1|) implies that the configurational entropy vanishes linearly in T at the 
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Kauzmann temperature Tk, and that Tk = Tq. The combination of Eqs. (j2.9J) and (|4.1|) 
yields 

AG*(T) ^ 7770*0 
fc B T ~3s c (T)' 1 • ] 

which is essentially the Adam-Gibbs formula near T . 

Definitions of the glass temperature T g generally have the form 

k B T g ^ 4 UJ v (T g — T ) " A " ^ 

where A 9 is a large number of order 30 or so, chosen roughly to represent the observable 
limits of long relaxation times or high viscosities. Then T g = To, and the fragility m |4JJ is 

m = -T— (*£2i) „ i m"" »_K MM 1 " (44) 



Thus the excitation-chain theory implies that glasses are fragile when Ti nt To, and/or 
when the disorder strength 70 is small. Returning to the thermodynamic formula, Eq. (j4.1j) . 
we find that the jump in the specific heat at T g = T is 

Here we recover the conjectured proportionality between Ac p and m, but with two material- 
specific parameters, so and 70, that might account for the observed scatter in the exper- 
imental data. According to its definition, so should scale with the "bead" number of the 
moleculeSjWhich ordinarily is factored out in obtaining the linear relation between Ac p and 
m. (See £ 9|-) Because it is basically a geometrical quantity, not involving energy scales, 
70 may be roughly a constant, of order unity, for most glassy materials. This analysis also 
implies that 

R*{T g ) « (-) ^. (4.6) 
Vtt/ 7o 

Thus the critical length scale R* at the glass temperature is predicted to be about 30 
molecular spacings, independent of the fragility. This prediction seems to be qualitatively 
consistent with results shown by Berthier et al. 47( ; but those authors report a substantially 
smaller length scale. 

Clearly, this thermodynamic analysis is incomplete and highly speculative. The 
excitation-chain theory and its proposed extension to thermodynamics both need to be 
explored further and tested experimentally. We also need to develop a theory of the 
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crossover between super-Arrhenius and Arrhenius behavior at Ta, and to understand 
how the present solidlike formulation crosses over to liquidlike mode-coupling theories at 
higher temperatures. 49, |3| And we have yet to address the important issue of stretched- 



exponential relaxation. (See jsij for an idea about how stretched exponentials might appear 
in STZ theory.) 

However, the deepest theoretical uncertainty is the meaning of Eq. (j4.1|) . If correct, this 
relation implies that ergodicity is broken not just below the Kauzmann temperature but, to 
a continuously increasing extent as T falls below Ta, throughout the super-Arrhenius region 
T < T < Ta- The language that I have used to support this conjecture is at best suggestive, 
and is far from being a systematic derivation of that equation. A list of unanswered questions 
makes the uncertainties abundantly clear. Is the domain structure to be taken literally, or 
is it just a way of talking about long-lasting correlations? Might there be some kind of long- 
range order - orientational or perhaps something even more subtle - inside the domains? 
Might the excitation chains occur predominantly in the more highly disordered boundary 
regions and, if so, might we need to return to the two-dimensional picture proposed in Q]? 
Is a "frozen" region permanently frozen? Presumably not, because the domain boundaries 
must diffuse at something like the Vogel-Fulcher relaxation rate. On the other hand, if all 
of the dynamically accessible states of the system consist of frozen domains separated by 
unfrozen boundaries, then Eq. (j4.1|) might be justified as the appropriate statistical average. 
How might such a picture be made into a quantitative, predictive theory? What theoretical 
tools might be useful? 
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